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Abstract 

On  Ice  Station  Weddell,  some  characteristics  of  the  neutral-stability  air-ice  drag 
coefficient  at  a  reference  height  of  10m  (Cbm  o) were  observed  that  had  not  been 
documented  before.  The  main  finding  was  that  wind-driven  snow  continually 
alters  the  sea  ice  surface;  the  resulting  snowdrifts  determine  how  large  Cdnio 
is.  In  particular,  this  report  describes  three  observations  and  attempts  to  explain 
them:  1)  CpNiois  near  1 .5  x  10-3  when  the  wind  is  well  aligned  with  the  drifted 
snow;  2)  Cdnio  is  near  2.5  x  10~3  when  the  wind  makes  a  large  angle  with  the 
dominant  orientation  of  the  snowdrifts;  3)  Cdni  o  can  increase  by  20%  if,  after 
being  well  aligned  with  the  drift  patterns,  the  mean  wind  direction  shifts  by  as 
little  as  20°.  To  investigate  this  behavior  of  Cdnio-  this  report  adapts  a  model 
developed  by  Raupach  that  partitions  the  total  surface  stress  into  contributions 
from  form  drag  and  skin  friction.  With  reasonable  choices  for  free  model 
parameters  and  with  littlefine-tuning,  this  physically  based  model  can  reproduce 
the  three  main  observations.  In  other  words,  the  model  seems  to  include  the 
basic  physics  of  air-ice  momentum  exchange.  This  modeling  implies  that  1 0- 
cm-high  sastrugi-like  roughness  elements,  rather  than  pressure  ridges,  sustain 
most  of  the  form  drag  over  compact  sea  ice  in  the  western  Weddell  Sea .  Lastly, 
the  report  speculates  on  what  the  observations  and  this  model  say  about  how 
to  parameterize  Cdnio  °ver  snow-covered  sea  ice. 
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A  Physically  Based  Model  of  the 
Form  Drag  Associated  with  Sastrugi 

EDGAR  L  ANDREAS 


INTRODUCTION 

The  neutral-stability,  10-m  air-ice  drag  coeffi¬ 
cients  (CDNio)  reported  by  Andreas  and  Claffey 
(1995;  henceforth,  AC95)  clearly  depend  on  how 
wind-driven  snow  alters  the  sea  ice  surface.  Fig¬ 
ure  1  shows  the  time  series  of  Ice  Station  Weddell 
(ISW)  Cdnio  values  and  hourly  averages  of  the 
wind  speed  and  direction  measured  at  5  m. 
These  CDnio  values  range  between  1 .27 x  10-3  and 
2.54  xlO-3  but  are  not  randomly  scattered  be¬ 
tween  these  limits.  Rather,  they  generally  cluster 
or  depend  coherently  on  time.  For  example,  in 
AC95  we  focused  on  the  events  on  Julian  days  64 
and  65  and  79  and  80.  In  Figure  1,  we  see  that,  on 
these  days,  the  5-m  wind  speed  was  typically 
8  m/s  or  higher,  the  wind  direction  was  fairly 
constant,  and  Cdnio  decreased  consistently  with 
time. 

Our  explanation  for  this  sequence  was  that 
winds  of  6-8  m/s,  and  higher,  redistribute  the 
mobile  snow:  snowdrifts  perpendicular  to  the 
wind  erode;  and  the  new  drifts,  built  parallel 
with  the  mean  wind,  streamline  the  surface.  Con¬ 
sequently,  if  these  eroding  winds  persist  and 
their  direction  remains  relatively  constant  for 
just  12  hours,  Cdnio  will  decrease  by  as  much  as 
30%.  If,  after  establishing  a  drift  pattern,  how¬ 
ever,  the  wind  shifts  by  as  little  as  20°,  it  will  face 
a  rougher  surface,  and  measured  Cdnio  values 
will  increase. 

There  are  several  other  events  in  Figure  1  that 
demonstrate  this  same  smoothing  in  the  pres¬ 
ence  of  a  moderate,  directionally  constant  wind 
(for  example,  see  days  86  and  87).  Thus,  the 
sequence  is  repeatable  and,  in  fact,  seems  to  be 
the  key  to  understanding  air-ice  momentum  ex¬ 
change  in  the  western  Weddell  Sea — and  proba¬ 
bly  over  most  other  sea  ice  regions  as  well. 

Form  drag  sustained  by  the  snowdrifts  must 
explain  the  behavior  of  the  Cdnio  values  that  we 
observed.  Others  have  also  recognized  that  form 
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Figure  1.  Time  series  of  the  Ice  Station  Weddell  (ISW) 
Qd^iq  values  reported  in  Andreas  and  Claffey  (1995) 
and  hourly  averages  of  the  wind  speed  and  direction  at 
5  m. 


drag  can,  at  times,  be  an  important  means  for 
transferring  momentum  to  a  surface.  For  exam¬ 
ple,  Lettau  (1969)  offered  a  heuristic  formula  in¬ 
volving  obstacle  geometry  and  silhouette  area  to 
parameterize  the  roughness  lengths  that  Kutz- 
bach  (1961)  measured  during  his  famous  bushel 
basket  experiment  on  Lake  Mendota.  Arya  (1973, 
1975)  assumed  that  pressure  ridges  were  the 
roughness  elements  responsible  for  most  of  the 
form  drag  over  Arctic  sea  ice  and  developed  a 
model  to  partition  the  total  stress  into  skin  fric¬ 
tion  and  form  drag  on  the  basis  of  ridging  inten¬ 
sity.  Banke  et  al.  (1976,  1980)  investigated  the 
pressure  distribution  around  Arctic  pressure 
ridges  1-3  m  high  in  an  attempt  to  determine  the 
form  drag  coefficient  for  these.  But  their  ultimate 
parameterization  for  Cdnio  (also  Banke  and 
Smith  1973)  is  more  in  line  with  our  findings; 
they  showed  that 

103  CDNio  =  1.10  + 0.072  §  (1) 

where  £  (in  centimeters)  is  the  root-mean-square 
surface  roughness  for  roughness  elements  with 
wavelengths  shorter  than  13  m.  That  is,  the 
smaller  scale  roughness  dominates  the  momen¬ 
tum  exchange.  Jackson  and  Carroll  (1978) 
reached  the  same  conclusion.  They  showed  that 
the  angle  between  the  mean  wind  and  the  domi¬ 
nant  axis  of  the  sastrugi  at  the  South  Pole  deter¬ 
mined  the  roughness  length;  they  then  devel¬ 
oped  a  model  based  on  Lettau's  to  explain  this 
effect. 

More  recently,  Raupach  (1992;  henceforth, 
R92;  also  Raupach  et  al.  1993)  developed  a  more 
sophisticated  model,  based  on  dimensional  anal¬ 
ysis  and  two  physically  reasonable  hypotheses, 
to  partition  the  total  surface  stress  into  contribu¬ 
tions  from  form  drag  and  from  the  stress  on  the 
underlying  surface.  That  model  is  the  basis  of  the 
analysis  here.  I  adapt  Raupach's  model  to  evalu¬ 
ate  the  stress  partitioning  for  a  field  of  sastrugi- 
like  snowdrifts.  That  model  predicts  the  drag  co¬ 
efficient  evaluated  at  the  top  of  the  roughness  el¬ 
ements;  I  show  how  this  quantity  is  related  to 
Cdnio/  which,  in  turn,  is  monotonically  related  to 
the  roughness  length  Zq  through 

Cdnio  =  fc2{In[(10-rf)/z0]}-2.  (2) 

Here,  k  (=  0.4)  is  the  von  Karman  constant,  and  Zo 
and  the  displacement  height  d  must  be  in  meters. 

My  objective  is  to  see  whether,  with  reason¬ 
able  choices  for  snowdrift  geometry  and  other 


free  model  parameters,  R92's  model  can  account 
for  the  behavior  of  Cdnio  at  ISW  that  we  docu¬ 
mented  in  AC95.  In  particular,  the  three  main  ob¬ 
servations  that  I  hope  to  explain  are  as  follows: 

1.  Cdnio  is  near  1.5  x  10-3  when  the  wind  is 
aligned  with  a  well-developed  field  of  drifted 
snow. 

2.  Cdnio  is  near  2.5  x  10-3  when  the  mean  wind 
direction  makes  a  large  angle  (say  50-140°)  with 
the  dominant  axis  of  the  drifts. 

3.  Cdnio  can  increase  by  roughly  20%  if,  after 
being  aligned  with  a  developed  (or  developing) 
field  of  drifted  snow,  the  wind  turns  by  as  little  as 
20°. 

RAUPACH'S  MODEL  ADAPTED  TO 
DRIFTED  SNOW 

Geometry 

Because  during  our  deployment  on  ISW  the 
wind  was  directionally  quite  variable,  the  sur¬ 
face  never  developed  the  deep,  elaborate  fields 
of  sastrugi  that  Mellor  (1965)  and  Jackson  and 
Carroll  (1978)  observed  on  the  Antarctic  conti¬ 
nent.  Rather,  on  ISW,  the  drifts  were  sparsely  dis¬ 
tributed  and  rudimentary.  When  the  wind 
reached  6-8  m/s,  linear  sastrugi-like  drifts 
would  begin  forming  behind  any  protruding 
obstacles.  The  drifts  would  grow  as  the  wind 
persisted;  but  when  the  wind  direction  changed, 
these  drifts  would  erode  in  winds  above  the  6-  to 
8-m/s  threshold.  Thus,  the  drifts  that  deter¬ 
mined  the  surface  roughness  on  ISW  were  gener¬ 
ally  small.  The  scour  marks  that  Allen  (1965)  de¬ 
scribed  and  named  scour-remnant  ridges  may 
well  be  precursors  to  these  rudimentary  sastrugi. 

Figure  2  shows  a  typical  leveling  profile  of  the 
snow  and  sea  ice  surfaces  in  the  vicinity  of  the 
profiling  mast  described  in  AC95.  We  made  sev¬ 
en  such  leveling  surveys  during  our  deployment 
on  ISW.  The  profiles  show  that  snow  tends  to  col¬ 
lect  around  obstacles  in  the  underlying  ice.  The 
root-mean-square  (rms)  roughness  of  the  snow 
surface  in  this  particular  profile  is  11  cm;  its  ice- 
surface  roughness  is  9  cm.  In  the  seven  surface 
profiles  that  we  collected,  the  rms  snow-surface 
roughness  ranged  from  9  to  26  cm;  the  ice-sur¬ 
face  roughness  ranged  from  9  to  28  cm.  A  value 
in  this  range  is  a  reasonable  choice  for  the  nomi¬ 
nal  height  of  the  snow-surface  roughness  ele¬ 
ments. 

Figure  3  shows  the  idealized  geometry  of  the 
roughness  elements  that  I  will  treat  with  R92's 
model.  These  are  rudimentary  sastrugi  with 
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Figure  2.  Snow  and  ice  surface  elevations  near  our  ISW  profile 
tower  on  Julian  day  78  (in  1992).  The  sampling  interval  was  0.5  m. 
The  profiling  line  is  not  necessarily  aligned  with  the  mean  wind. 
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Figure  3.  Geometry  of  the  sastrugi-like  snowdrifts  being  modeled  here. 


height  h,  width  mh,  and  length  nh.  The  angle  <1> 
describes  how  the  mean  wind  is  oriented  with 
respect  to  the  sastrugi.  When  the  wind  has  been 
blowing  with  constant  direction  at  speeds  of 
6-8  m/s,  or  higher,  for  several  hours,  <E>  will  be 
zero. 

From  the  discussion  in  the  penultimate  para¬ 
graph,  we  see  that  h  is  5-20  cm.  Descriptions  of 
drifted  snow  in  the  literature  (Finney  1939,  Mel- 
lor  1965,  Jackson  and  Carroll  1978)  and  Plate's 
(1971,  Figure  4.17)  pressure  distribution  around 
roughness  elements  similar  to  that  in  Figure  3 
suggest  that  n  is  about  10.  From  Mellor's  (1965) 
description  of  sastrugi,  I  assume  that  m  is  about 
4.  Notice,  these  parameters  imply  that,  in  Figure 
3 


P  =  arctan(m/2n)  =  11.3°.  (3) 

One  of  the  key  parameters  in  R92's  model  is  X, 
the  frontal  area  of  a  typical  roughness  element 


per  unit  ground  area.  Suppose  the  fraction  of  the 
area  covered  by  the  sastrugi-like  roughness  ele¬ 
ments  is  y.  Thus,  from  Figure  3,  the  unit  area  repre¬ 
sented  by  each  roughness  element  is  mnh2/ 2y.  The 
base  length  of  the  triangular  frontal  area  that  each 
roughness  element  presents  to  the  wind,  L(<I>),  is  a 
function  of  <£>.  Likewise,  the  silhouette  area  pre¬ 
sented,  A(<h),  is  simply 


A(o)=hm. 

Consequently,  X  is  also  a  function  of  O 
mnh L  mnh 


(4) 

(5) 


In  Figure  3,  we  see  that  there  are  four  O  regions 
for  which  a  roughness  element  presents  different 
faces  to  the  wind  for  0  <  O  <  180°  (symmetric  for 
0  >  <!>  >  -180°).  These  and  the  associated  values  of 
L,  A  and  X  follow. 
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Region  I;  0°  <  3>  <  (3 

Here  the  wind  sees  only  the  front  triangular 
face.  Thus 


L(®)  =  mh  cos® 

(6a) 

A(®)  =  (mh2/2)  cos® 

(6b) 

X(®)  =  (y/n)  cos®  . 

(6c) 

Region  II:  (3  <  O  <  90° 

Here  the  wind  sees  both  the  front  triangular 
face  and  the  side  ridge.  Thus 

L(®)  =  nh  sin®  +  (mh/ 2)  cos® 

(7a) 

A(®)  =  (h2/2)[n  sin®  +  (m/2)  cos®] 

(7b) 

X(®)  =  y[m_1  sin®  +  (2 n)-1  cos®] . 

(7c) 

Region  III:  90°  <  ®  <  180°-p 

Here  the  wind  sees  only  the  side  ridge.  Thus 

L(0)  =  nh  sinO  -  (mh/2)  cosO 

(8a) 

A(O)  =  (h2/2)  [n  sinO  -  (m/2)  cosO] 

(8b) 

A,(0)  =  yfw1  sinO  -  (2 n)-1  cosO]  . 

(8c) 

Region  IV:  180°  -  (3  <  O  <  180° 

Here  the  flow  sees  both  side  ridges 
approaches  from  the  rear  of  the  drift.  Thus 

as  it 

L(<1>)  =  -mh  cos<2> 

(9a) 

A(<D)  =  -(mh2/2)  cosO 

(9b) 

X(<&)  =  -(y/n)  cos<l> . 

(9c) 

Partitioning  the  stress 

R92's  intent  was  to  partition  the  total  surface 
stress  (x)  into  contributions  from  form  drag  (xr) 
and  from  stress  on  the  underlying  surface  (the 
skin  friction,  Xg).  In  general 

x  =  xR+xs.  (10) 

There  is  another  potential  sink  for  the  momen¬ 
tum  in  eq  10  that  R92  did  not  treat.  Over  erodible 
surfaces,  saltating  particles  may  extract  momen¬ 
tum  from  the  air  because  of  their  inertia.  Neither 
xR  nor  Xg  reflects  this  momentum  exchange.  That 
saltation  layer,  however,  is  quite  thin — on  the 
order  of  a  centimeter  (Owen  1964,  Radok  1968, 


AC95).  Thus,  since  R92  models  eq  10  for  heights  of 
order  h  and  higher — well  above  the  saltation 
layer — eq  10  need  not  include  a  term  to  account 
for  particle  inertia.  Above  the  saltation  layer,  the 
concentration  of  blowing  and  drifting  snow  is  too 
low  to  affect  the  dynamics  of  the  flow. 

R92  modeled  the  stress  on  the  underlying  sur¬ 
face  as 


Ts(^)  =  pCShL^exp 


-cxA,(®)^ 

u 


(11) 


Here,  Cg^  is  the  drag  coefficient  of  the  underlying 
surface  referenced  to  a  height  h,  which  is  where 
the  wind  speed  iih  is  evaluated.  Also  in  eq  11,  p  is 
the  air  density;  u the  friction  velocity,  w3f=(x/p)1/2; 
and  C\,  an  empirical  constant.  The  exponential 
term  in  eq  11  accounts  for  the  sheltering  of  the 
underlying  surface  by  the  roughness  elements. 
That  is,  when  the  roughness  elements  are  small  or 
sparse  and  X(O),  thus,  is  small,  eq  11  reduces  to 
xs  =  pCshLZh,  the  usual  expression  for  the  skin  fric¬ 
tion  over  a  smooth  surface.  On  the  other  hand, 
when  the  roughness  elements  have  a  large  frontal 
area  or  are  densely  packed,  X(O)  is  large,  and  xs 
approaches  zero.  In  this  case,  the  roughness  ele¬ 
ments  completely  shelter  the  surface;  thus,  skin 
friction  can  provide  none  of  the  stress. 

We  can  infer  the  value  of  CSh  appropriate  for 
smooth,  snow-covered  surfaces  from  Overland's 
(1985)  review  and  from  the  measurements  report¬ 
ed  by  Banke  et  al.  (1980)  and  Kondo  and  Yamaza- 
wa  (1986).  The  lowest  value  for  Cdnio  reported  by 
Overland  and  Kondo  and  Yamazawa  is  roughly 
1.1  x  10-3.  Coincidentally,  Figure  4  in  Banke  et  al. 
implies  that  the  value  of  Qdnio  for  completely 
smooth  sea  ice  is  1.10  x  10-3.  Call  this  minimum 
skin  friction  coefficient  CgiQ. 

Converting  Cgio  to  Cgh,  however,  is  not  straight¬ 
forward:  We  do  not  know  the  displacement  height 
d  a  priori;  and  near  the  tops  of  the  roughness  ele¬ 
ments,  the  semi-logarithmic  profile  law  (see  eq  2) 
is  not  accurate.  I  will  explain  how  I  deal  with  these 
complications  later;  here,  suffice  it  to  say,  I  obtain 
Csh  from  CS10. 

R92  wrote  a  general  expression  for  the  force 
caused  by  form  drag  on  an  isolated  roughness  ele¬ 
ment  as 


FD=pCRAUl  (12) 

where  Cr  is  the  drag  coefficient  of  the  obstacle, 
and  A  is  its  silhouette  area.  Here,  I  must  modify 
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Raupach's  model,  however.  In  my  model,  the  ob¬ 
stacle's  aerodynamic  properties  likely  change  as 
its  orientation  with  the  wind  changes.  Other¬ 
wise,  the  concept  of  streamlining  would  be 
meaningless.  Thus,  Cr  and  A  must  depend  on  O. 
But  because  I  could  find  no  study  in  the  literature 
that  reported  the  aerodynamic  properties  of  an 
obstacle  like  that  depicted  in  Figure  3, 1  assume 
that  the  drag  contributions  from  the  various  fac¬ 
es  are  additive.  Mathematically,  my  hypothesis  is 

FD(0)  =  pUh2XCRi(<I»)Ai(4»)  (13) 

i 

where  Cr*  is  the  drag  coefficient  for  a  particular 
face  of  the  sastrugi-like  obstacle,  and  A[  is  the  sil¬ 
houette  area  of  that  face.  The  sum  is  over  all  faces 
presented  to  the  wind  at  angle  O. 

The  obstacle  in  Figure  3  presents  only  three 
faces.  When  0°  <  O  <  90°,  the  wind  can  see  the 
front  triangular  face;  this  has  form  drag  coeffi¬ 
cient  Cri.  When  p  <  <£  <  180°  -  P,  the  wind  can  see 
the  side  ridge;  this  has  drag  coefficient  Cr2. 
When  180°  -  p  <  O  <  180°,  the  wind  can  see  only 
the  rear  silhouette;  this  has  drag  coefficient  Cr3. 

As  I  mentioned,  I  found  no  study  that  evaluat¬ 
ed  CRfor  the  obstacle  geometry  I  use  here.  But  on 
reviewing  several  studies  that  treated  somewhat 
similar  geometries  (i.e.,  Arie  and  Rouse  1956; 
Arya  1973,  1975;  Banke  et  al.  1976,  1980;  Taylor 
1988),  I  estimated  what  I  feel  are  representative 
values  for  the  three  faces.  Here,  I  use  Cri  =  0.10, 
Cr2  =  0.30,  and  Cr3  =  0.30. 

Continuing  now  with  R92's  model,  if  p  rough¬ 
ness  elements  cover  a  ground  area  S,  the  stress 
they  produce  through  form  drag  is 


p 


vY 

1_SfcJ  •  (14> 


Here,  V  is  the  volume  sheltered  by  an  individual 
roughness  element;  the  term  raised  to  power  p, 
thus,  represents  a  mutual  sheltering  effect. 

In  my  model,  from  eq  5 


V  _  2y  _  1(0) 

S  mnh2  A(O)  ’ 


(15) 


Using  this  in  eq  14  gives 


m 

A(d>) 


p^hlQ^o)^®) 

i 


i(®)y 

phA(0)  J 


W 


(16) 


In  the  limit  of  large  p 

-X(<Z>)V 
hA(<t>) 


!  mv 

phA{<t>)  j 


-exp 


(17) 


Finally,  following  R92,  for  m  and  n  greater  than  1, 
I  estimate  the  sheltered  volume  to  be 


V  =  ±c2L(4>)h2^  (18) 

z  u* 

where  c2  is  another  empirical  constant.  Hence, 
using  eq  5, 17  and  18  in  eq  16,  we  get 


xR(®)  =  1(0)1 


puffier, 


exp 


_c2  1(®)^|. 


(19) 


In  eq  19,  define 


(20) 

For  the  four  O  regions  identified  earlier,  CR  fol¬ 
lows. 


Region  I:  0°  <  O  <  p 

Cri  =  CR1.  (21a) 

Region  II:  P  <  ®  <  90° 


£  _  CR1m  cos  O  +  Cm  [n  sin  Q-(m/ 2)  cos  O] 

1511  nsin®  +  (m/2)cos® 

(21b) 

Region  III:  90°  <  O  <  180°  -  P 

Qun  =  0*2  •  (21c) 

Region  IV:  180°  -  p  <  $  <  180° 

Criv  =  CR3-  (21d) 

Combining  eq  11, 19  and  20  in  eq  10,  we  get 


Q_ 

II 

sP'r 

CSh  exp 

si 

g 

O 

1 

L 

+l(0)CR(0)exp 


-c2l(0) 


Uh 


(22) 
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To  make  solving  eq  22  easier,  R92  assumed  that 
Cj  =  c2  =  c;  for  c  he  used  0.25,  which  provided  the 
best  fit  to  the  data  that  he  tested  his  model 
against.  With  this  assumption,  eq  22  reduces  to 


x(O)  =  plfh[cSh  +X(0)CR(0)]exp 


-cX(O) 


Uh 


(23) 


Notice 


Displacement  height 

To  convert  the  total  drag  coefficient  CDh  to  its 
10-m  value,  we  need  to  know  the  wind  speed 
profile  law.  When  a  surface  is  covered  with 
roughness  elements  that  are  much  larger  than  z0, 
it  is  often  necessary  to  include  a  displacement 
height  d  (Thom  1971)  in  the  semi-logarithmic 
wind  profile.  That  is 


t(<I>)_m*  =J_ 

,  r7  ,t2  ~  ”  2 

p^h  uh  r 


(24) 


is  the  drag  coefficient  for  the  total  surface  stress 
referenced  to  height  h.  On  inserting  this  defini¬ 
tion  and  rearranging  terms  in  eq  23,  we  get 


cl(0)r 


exp 


-cX(0)r 

2 


:^M[cSh  +  X(0)CR(0)]'1/2.  (25) 


This  has  the  form 
X  e~x  =  a 


(26) 


u*  k  l  zo 


(29) 


where  li(z)  is  the  average  wind  speed  at  height  z. 

But  near  the  tops  of  the  roughness  elements, 
eq  29  cannot  be  accurate  because  of  the  enhanced 
turbulence  there  created  by  flow  around  these 
obstacles.  R92  therefore  hypothesized  that,  for 
h<z<zw 


1 

k 


f 

In 


z-d 


f 

+v 

V 


z-d  1 

Zv-dJ 


(30) 


where  \ji  is  a  profile  influence  function  similar  to 
the  stability  corrections  required  by  Monin- 
Obukhov  similarity.  At  zw  and  above,  \|/  is  zero. 
R92  went  on  to  show  that 


where  a  is  a  constant  for  a  given  value  of  <t>. 

I  solved  eq  26  for  X  using  Newton's  method; 
the  solution  required  only  two  to  three  iterations. 
There  are  multiple  solutions  to  eq  26,  however; 
we  must  select  the  right  one.  As  R92  explained, 
eq  26  has  the  following  properties: 

For  0  <  a  <  e-1,  eq  26  has  two  solutions,  Xi  and 
X2,  where  Xj  <  1  and  X2  >  1. 

For  a  =  e-1,  X  =  1. 

For  a  >  e-1,  there  is  no  solution  to  eq  26. 

The  only  nontrivial  case  is  the  first  of  these;  and 
Xi  is  the  only  physical  solution,  since  X  must  ap¬ 
proach  zero  as  either  A,(0)  or  a  approaches  zero. 
Thus,  to  start  Newton's  method,  my  first  guess 
was  X  =  a. 

Solving  eq  26  yields  X,  which,  in  turn,  is  relat¬ 
ed  to  the  total  drag  coefficient  we  seek 


r- 1/2  _  2X 

Dh  a(4>)‘ 


(27) 


Notice,  to  check  this  solution,  we  can  make  an 
asymptotic  approximation  in  eq  23.  When  A,(<£>)  is 
near  zero 


CDh=CSh  +  X(0)CR(0).  (28) 


V 


f  z-d' 


VZw 


Zw-d 


and  that 


d  —  6?) 

where  cw  is  a  constant.  He  also  defined 


Vh  =  V 


(  h-d  ' 

VZW  ~~  d  j 


—  ln(cw )  +  cw  1 


(31) 


(32) 


(33) 


which  will  be  useful  shortly.  Notice,  there  are  er¬ 
rors  in  two  of  R92's  equations  (eq  29  and  31).  In 
an  unpublished  correction  to  his  paper  dated  23 
October  1992,  Raupach  explained  that,  in  light  of 
these  errors,  cw  is  approximately  4  rather  than 
1.5,  the  value  he  originally  gave.  Hence, 
\ph  =  0.64,  a  value  independent  of  surface  proper¬ 
ties. 

Equation  32  implies  that 

zw  =cwh  +  (l-cw)d.  (34) 

With  cw  =  4  and  the  assumption  that  zw  >  h,  eq  34 
implies  that  d  <  h:  The  displacement  height  is  al¬ 
ways  less  than  the  height  of  the  roughness  ele¬ 
ments. 
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To  obtain  Cdnio  from  Cph/ we  must  match  the 
profile  laws,  eq  29  and  30,  at  zw.  From  eq  29  and 
32,  we  have  the  following  results  for  z  >  zw 


r-i/2  _U(  10) 

CDN10  =  ~ 


10  -P 

zo  , 


and 


(35) 


Notice,  from  eq  19, 20  and  23 

%_  MofoW  (44) 

t  Csh+X(<S)CR(4>) 

R92  then  related  dR  to  h .  Following  his  argu¬ 
ments  but  applying  them  to  the  specific  geome¬ 
try  of  my  model,  I  get 


C-l/2  =ii(Zw)  1  ^ 


+  ln(cw)  (36) 


where  h,  d,  zq  and  zw  must  be  in  meters.  Subtract¬ 
ing  eq  35  from  36,  we  get 


UDzw  “UDN10 


From  eq  30-33,  for  h  <  z  <  zw 


l  +  ¥hl  (38) 

w*  Kl  \  zo  J 

and 

<39> 

U *  KL  V  Z0  ) 

Subtracting  eq  38  from  39  yields 

C^=C£jh/2+i[ln(cw)-¥h].  (40) 

On  equating  eq  37  and  40,  we  finally  can  relate 
Cdnio  to  Coh 


dR=h-cdC^2[A(0)/X(0)f/2  (45) 

in  which  c d  is  a  constant  that  R92  took  to  be  0.6. 
Inserting  eq  5  in  this  and  combining  it  with  eq  43 
yields 

d  =  h(zR/  t)[l  -  cjCp/2  ( mn  /2y)1/2].  (46) 

Thus,  it  is  easy  to  find  Tr/t  from  eq  44;  and  once 
my  model  yields  Cp h,  it  is  simple  to  compute  d. 
Using  this  value  in  eq  41  finally  gives  the  quan¬ 
tity  we  seek,  Cdnio- 

I  pointed  out  in  the  previous  section  that  we 
need  to  convert  Csio  to  Csh  before  we  can  begin 
computations.  Equation  41  also  makes  this  con¬ 
version.  Unfortunately,  we  need  Csh  before  we 
have  obtained  d .  I  could  handle  this  problem  by 
iterating  the  entire  set  of  equations  on  d.  But,  as  I 
will  show  in  the  next  section,  d/h  *  0.3,  a  result 
also  consistent  with  eq  46.  At  this  step  in  the 
computations,  this  simple  approximation  is  rea¬ 
sonable — especially  in  light  of  uncertainties  in 
the  other  model  parameters — because  Cdnio  is 
not  very  sensitive  to  Csio  for  the  range  of  values 
that  this  parameter  can  realistically  assume. 


C™10  -  CDh2  +  \  )  Vh  •  {41) 

We  still  do  not  know  the  displacement  height  d, 
however. 

Thom  (1971)  identified  d  as  the  effective  level 
at  which  the  roughness  elements  absorb  the  mo¬ 
mentum  being  transferred  to  the  surface.  R92 
used  this  definition  to  derive  d.  That  is,  if  ^r  is  the 
centroid  of  the  form  drag  on  the  roughness  ele¬ 
ments  and  if  dg  is  the  centroid  at  which  the  skin 
friction  acts 

dx  =  dRxR  +  dgXg  z11  dRxR  (^2) 

because  dg  =  0  by  definition.  Thus 

d  =  dR(xR/t).  (43) 


RESULTS 

I  have  not  yet  discussed  the  value  of  y,  the 
fractional  coverage  of  sastrugi-like  roughness  el¬ 
ements.  Vladimir  Churun*  tried  to  quantify  the 
roughness  of  the  ISW  floe  using  the  radar  on  the 
Akademik  Fedorov  at  the  time  the  station  was  de¬ 
ployed.  His  survey  suggested  that  hummock 
coverage  varied  over  the  floe  from  10  to  30%. 
While  hummock  coverage  is  not  the  same  as  sas- 
trugi  coverage,  my  personal  experience  on  ISW 
suggests  that  10-30%  sastrugi  coverage  is  also 
about  right.  Notice,  because  of  the  geometry  of 
the  sastrugi  that  I  am  modeling  (Fig.  3),  the  tight¬ 
est  possible  packing  of  roughness  elements  will 
yield  a  y  value  of  only  0.5.  Thus,  my  first  guess, 

*  Personal  communication  (Arctic  and  Antarctic  Research 
Institute,  St.  Petersburg,  Russia). 


y  =  0.15,  represents  a  conservative  but  not  insig¬ 
nificantly  small  sastrugi  coverage. 

Figure  4  shows  two  sets  of  predictions  for 
Cdnio  as  a  function  of  based  on  the  model 
developed  in  the  preceding  section.  Figure  4a 
shows  computations  of  Cdnio  for  y  =  0.15  and  for 
sastrugi  heights  of  5, 10  and  20  cm,  typical  heights 
represented  in  our  snow-surface  roughness  sur¬ 
veys  (e.g..  Fig.  2).  The  10-cm  case  fits  our  observa¬ 
tions  quite  well.  Here,  Cdnio  is  1.43-1.45  xlO-3 
when  the  flow  is  head-on  or  within  ±12°  of  head- 
on  (that  is,  when  <E>  <  (J).  These  values  agree  well 
with  the  nominal  minimum  value  for  Cdnio  of 
1.5  x  10~3  that  we  reported  in  AC95. 

When  <X>  gets  larger  than  11.3°  (see  eq  3),  the 
wind  can  see  the  side  ridge  of  the  obstacle.  In  ef¬ 
fect,  it  sees  a  larger  silhouette  area  and  a  higher 
value  of  Cr;  Cdnio  roust  thus  increase.  Figure  4 
shows  that  if  the  wind  rotates  just  20°  from  head- 
on,  Cdnio  will  increase  by  15-24%,  in  good  agree¬ 
ment  with  our  observations  in  AC95. 

As  <E>  increases  and  the  wind  thus  has  an  in¬ 


creasing  component  perpendicular  to  the  axis  of 
the  sastrugi,  Cdnio  continues  increasing.  Cdnio  bas 
its  maximum  at  O  =  101.3°,  where  it  is  2.73  x  10~3 
for  the  y  =  0.15  and  h  =  10  cm  case.  Again,  m  and  n 
set  the  location  for  this  maximum;  this  is  where 
A($)  has  its  maximum.  The  Cdnio  value  there  is  in 
fair  agreement  with  the  maximum  value  that  we 
observed  on  ISW,  2.54  x  10~3  (AC95). 

As  <D  increases  toward  180°,  the  model  is  entire¬ 
ly  speculative,  since  we  have  few  data  with  which 
to  compare  such  a  reversal  in  the  wind  and  be¬ 
cause  Cr3  is  the  least  certain  of  the  form  drag  coef¬ 
ficients.  For  now,  the  decrease  as  O  approaches 
180°  seems  reasonable  because  the  silhouette  area 
is  decreasing  rapidly.  But  the  absolute  values  for 
Cdnio  hr  this  region  are  certainly  open  to  debate. 

In  Figure  4a,  the  calculations  with  h  values  of  5 
and  20  cm  are  less  realistic  than  for  h  =  10  cm.  For 
h  =  5  cm,  the  maximum  Cdnio  value  calculated  is 
8%  less  than  our  observed  maximum;  for  h  =  20 
cm,  the  calculated  maximum  is  27%  higher  than 
what  we  observed.  For  all  three  heights,  however. 


Figure  4.  Model  calculation  of  Cdnio  as  a  function  of  the  wind  orientation  (O)  with  respect  to  the 
dominant  axis  of  the  sastrugi.  Model  parameters  are  m  =  4,  n  =  10,  C sio  =  1-10  x  10“3,  CR1  =  0.10, 
Cr2  =  0.30,  Cr3  =  0.30,  c  =  0.25 ,  cw  =  4  and  c ^  =  0.6. 
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Cdnio  values  for  head-on  flow  are  in  a  reasonable 
range,  1.36-1.51  x  10‘3. 

Figure  4b  investigates  the  sensitivity  of  the 
model  to  y  when  ft  =  10  cm.  For  this  value  of  ft, 
y  =  0.10  yields  Cdnio  values  for  head-on  flow  of 
1.36  x  10-3.  Although  we  did  measure  four  Cdnio 
values  in  this  range,  the  more  likely  lower  limit  is 
about  1.5  x  10~3.  Also,  as  O  increases,  the  predict¬ 
ed  Cdnio  values  reach  a  maximum  that  is  6%  be¬ 
low  our  observed  maximum.  When  y  is  0.20  and 
0.25,  the  predicted  head-on  values  for  Cdnio  are 
roughly  correct,  1.50-1.55  x  10-3.  But  for  these  y 
values,  the  maximum  predicted  values  are,  re¬ 
spectively,  17  and  26%  above  our  observed  maxi¬ 
mum.  Notice,  y  =  0.25  means  that  the  sastrugi 
cover  half  the  area  that  they  would  if  they  were 
packed  as  tightly  as  possible.  This  value  is,  there¬ 
fore,  probably  an  unrealistically  large  sastrugi 
coverage  for  ISW. 

For  the  y  =  0.15  and  ft  =  10  cm  case — the  one 
that  fits  the  data  best — I  also  ran  the  model  with 
Csio  =  1.00  x  10-3  (not  shown),  roughly  the  small¬ 
est  value  compatible  with  the  surveys  I  cited  ear¬ 
lier.  Resulting  predictions  for  Cdnio  were  4-6% 
lower  for  0°  <  O  <  20°  and  for  160°  <  O  <  180°  than 
for  calculations  based  on  Csio  =  1.10  x  10-3.  For 


Figure  5.  Stress  partitioning ,  Tr/t,  versus  the  mean 
wind  direction  with  respect  to  the  dominant  axis  of 
the  sastrugi.  Realize  that  %$/%  =  1-tr/ t.  Here ,  y  is  the 
fractional  sastrugi  coverage ,  and  h  is  the  sastrugi  height 
in  centimeters. 


Figure  6.  Computed  displacement  height  d  scaled  by 
sastrugi  height  h  versus  the  direction  of  the  mean 
wind  with  respect  to  the  sastrugi.  Here ,  y  is  the  frac¬ 
tional  sastrugi  coverage ,  and  h  is  in  centimeters. 


30°  <  <X>  <  160°,  the  differences  were  only  2-3%. 
Consequently,  the  choice  of  Csio  does  not  have  a 
big  effect  on  the  results,  as  long  as  the  value  is  a 
reasonable  one.  Secondly,  my  handling  of  the  dis¬ 
placement  height  in  converting  from  Csio  to  Csh  is 
appropriate,  given  this  weak  sensitivity  to  Csio- 

Figure  5  shows  some  examples  of  how  the 
stress  is  partitioned  between  form  drag  and  skin 
friction.  I  plot  Tr/ t;  clearly,  tgA^l  “  tr/ ^  the 
figure,  the  (y  =  0.10,  h  =  10  cm)  and  the  (y  =  0.20, 
h  -  20  cm)  lines,  respectively,  show  the  minimum 
and  maximum  values  of  Tr/t  that  I  obtained  for 
the  range  of  model  parameters  considered.  My 
computations  suggest  that,  even  for  a  very  rough 
surface,  y  =  0.20  and  h  =  20  cm,  skin  friction  pro¬ 
duces  more  than  half  of  the  total  stress  when  the 
flow  is  approximately  head-on.  In  other  words, 
the  streamlined  shape  of  the  sastrugi  is  quite  ef¬ 
fective  at  reducing  form  drag.  When  the  flow  is  at 
right  angles  to  the  sastrugi,  however,  form  drag 
dominates  the  total  stress  because  the  wind  has 
so  much  more  surface  to  push  against.  This  dom¬ 
inance  of  form  drag  for  such  wind  orientations 
explains  why  the  drifts  erode  so  quickly  (AC95). 

Equation  46  yields  the  displacement  height; 
Figure  6  shows  several  computations  of  d/ft .  For 
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densely  packed  roughness  elements,  such  as 
crops,  trees  or  other  vegetation,  d/his  typically 
reported  to  be  about  0.7  (Stanhill  1969,  Thom 
1971,  Fazu  and  Schwerdtfeger  1989).  For  the 
more  widely  spaced  sastrugi  that  I  am  consider¬ 
ing,  however,  d/h  is  0.15-0.35.  This  is  also  rough¬ 
ly  what  R92  (his  Figure  5)  computed  for  X  values 
in  the  range  relevant  to  this  problem,  0.01-0.06. 

These  plots  justify  my  use  of  d/h  =  0.3  when  I 
had  to  convert  Csio  to  C$ h  in  the  section  on  Parti¬ 
tioning  the  Stress .  They  also  show  that,  for  sastru¬ 
gi  typically  10  cm  high,  d  is  usually  less  than  3 
cm.  We  were  thus  safe  to  ignore  d  in  analyzing 
our  wind  speed  profiles  in  AC95.  For  the  lowest 
profiling  level  used  in  AC95,  0.5  m,  including  a 
displacement  height  of  3  cm  in  a  computation 
based  on  eq  29  would  yield  a  u*  value  roughly 
1%  larger  than  a  computation  with  d  =  0.  For  the 
higher  levels,  the  difference  would  be  even  less. 
Since  the  nominal  uncertainty  in  our  wind  speed 
measurements  was  ±5  cm/ s,  this  small  bias  in 
our  computations  is  negligible. 

Figure  6  shows  that  d  has  some  interesting 
structure — it  is  not  a  constant  for  all  wind  direc¬ 
tions.  Equation  46  explains  this  behavior:  d/h  re¬ 
sponds  mainly  to  two  variables,  Tr/ x  and  Cd h; 
d/h  increases  linearly  with  Tr/t  but  decreases 
linearly  with  Cq^2.  The  peak  in  d/h  near  20-30° 
corresponds  to  the  rapid  increase  in  xR/ x  in  this 
range  (Fig.  5),  while  the  broad  minimum 
centered  near  100°  is  where  CD h  has  its  maxi¬ 
mum  (Fig.  4).  Another  curious  feature  in  Figure  6 
is  that  the  d  values  in  this  broad  minimum  are 
usually  less  than  the  values  for  head-on  flow. 
Clearly,  d  does  not  change  in  concert  with  z0  or 
CDh — parameters  that  quantify  the  aerodynamic 
roughness  of  the  surface.  This  is  an  important 
conclusion  because  others  have  claimed  that  d  is 
a  linear  function  of  zq  (e.g.,  Brutsaert  1982,  p.  113; 
Sugita  and  Brutsaert  1990). 

In  summary,  the  Raupach's  Model  Adapted  to 
Drifting  Snow  section  explains  that  my  adapta¬ 
tion  of  R92's  model  has  11  adjustable  parameters: 
h,  m,  n,  y,  CsiO/  Cri,  Cr2,  Cr3,  c,  cw  and  c With 
reasonable  choices  for  these  and  with  little  fine- 
tuning,  I  have  demonstrated  that  that  model  fits 
our  ISW  observations  quite  well.  I  do  not  mean 
to  imply  here,  however,  that  the  model  is  infalli¬ 
bly  accurate — only  that  it  seems  to  capture  the 
essential  physics  of  how  wind  transfers  momen¬ 
tum  to  snow-covered  sea  ice. 


DISCUSSION 

Primary  roughness  elements 

The  obvious  (and,  perhaps,  surprising)  conclu¬ 
sion  from  the  calculations  described  in  the  last 
section  is  that  roughness  elements  with  heights  of 
only  10  cm  can  have  a  dominant  effect  on  Cdnio* 
This  is  essentially  the  same  conclusion  that  Jack- 
son  and  Carroll  (1978)  reached.  Arya's  (1973, 
1975)  model,  in  contrast,  predicted  that  pressure 
ridges  1  to  4  m  high  accounted  for  most  of  the 
form  drag  that  determined  the  large-scale  rough¬ 
ness  and  drag  coefficient  over  Arctic  sea  ice. 

My  calculations  provide  a  theoretical  basis  for 
the  empirical  result  that  Banke  et  al.  (1976, 1980) 
obtained,  eq  1.  Their  roughness  parameter  £,  re¬ 
flects  the  small-scale  roughness;  they  reported  £ 
values  of  3-13  cm,  where  ^  came  from  integrating 
the  roughness  spectrum  over  all  wavelengths 
shorter  than  13  m  (also  see  Andreas  et  al.  1993). 
Such  wavelengths  are  far  smaller  than  either  Arc¬ 
tic  or  Antarctic  ridge  spacings  (Lytle  and  Ackley 
1991).  Thus,  eq  1  and  my  calculations  point  to  the 
same  conclusion:  Pressure  ridges  are  relatively 
unimportant  in  setting  the  local  drag  coefficient; 
Cdnio  responds  primarily  to  sastrugi-size  rough¬ 
ness  elements.  Pressure  ridges,  however,  will 
probably  be  important  in  establishing  an  areally 
averaged  or  effective  roughness  length,  Z0.  Fied¬ 
ler  and  Panofsky  (1972),  Arya  (1975),  Overland 
(1985),  and  Claussen  (1991),  among  others,  have 
offered  some  thoughts  on  inferring  Zq  from  z0. 
But  pursuing  that  connection  is  beyond  my  scope 
here. 

Parameterizing  the  drag  coefficient 

These  model  results  and  the  behavior  of  Cdnio 
over  snow-covered  sea  ice  that  we  documented  in 
AC95  suggest  what  might  be  necessary  in  a 
scheme  for  parameterizing  Cdnio-  The  snow  sur¬ 
face  is  not  in  general  isotropic,  especially  during, 
or  for  some  time  after,  high  wind  events.  The 
wind  builds  sastrugi  and  therefore  streamlines 
the  surface  in  the  mean  wind  direction.  Figure  4 
shows  how  Cdnio  behaves  for  various  wind  di¬ 
rections  once  the  surface  has  been  streamlined. 

Thus,  the  key  is  to  estimate  the  dominant  di¬ 
rection  of  the  sastrugi  without  the  benefit  of  in 
situ  observations.  Maybe  in  the  future,  satellites 
will  provide  real-time  information  on  ice-surface 
roughness,  including  any  preferential  alignment 
of  the  roughness  elements.  But  for  now,  all  we  can 
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do  is  estimate  the  orientation  of  roughness  ele¬ 
ments  by  tracking  the  history  of  surface  winds  that 
have  been  inferred  from  analyzed  geostrophic 
winds.  According  to  AC95,  if  the  surface  wind 
blows  at  8  m/s  or  higher  with  constant  direction,  it 
will  build  sastrugi  parallel  to  it  and,  thus,  stream¬ 
line  the  surface.  If  such  conditions  persist  for  12 
hours,  Cdnio  will  fall  from  its  initial  value  to  about 
1.5  x  10 ~3  and  will  remain  here  as  long  as  the  wind 
remains  aligned  with  the  sastrugi. 

If  subsequent  surface  winds  fall  below  6  m/s. 
Figure  4  shows  how  Cdnio  will  vary  as  the  wind 
turns  with  respect  to  the  axis  of  the  sastrugi. 

The  difficult  part  of  the  parameterization  will 
be  treating  directionally  variable  winds  of  6-8nvfe, 
or  higher.  These  are  strong  enough  to  begin  re¬ 
working  the  surface — eroding  the  snowdrifts 
present — but,  because  of  their  variability,  will  not 
necessarily  build  new  drifts  and  thereby  stream¬ 
line  the  surface  in  another  direction.  Neither  the 
data  in  AC95  nor  the  model  developed  here  pro¬ 
vides  a  clear  suggestion  on  how  to  evaluate  Cdnio 
for  such  winds. 

In  summary,  except  for  the  complication  just 
mentioned,  it  seems  possible  to  now  estimate 
Cdnio  over  snow-covered  sea  ice  by  tracking  the 
history  of  the  surface  wind  computed  from  ana¬ 
lyzed  pressure  fields.  I  have  described  the  algo¬ 
rithm  in  words;  coding  and  testing  it  is  the  next 
step. 

CONCLUSIONS 

I  have  adapted  R92's  model  for  the  form  drag 
on  a  field  of  interacting  roughness  elements  and 
applied  it  to  rudimentary  sastrugi-like  snow¬ 
drifts.  The  key  step  was  treating  the  complex  ge¬ 
ometry  of  the  sastrugi.  In  the  absence  of  laborato¬ 
ry  values  for  the  form  drag  on  the  various  faces  of 
the  sastrugi,  I  assumed  that  the  form  drag  reflects 
additive  contributions  from  the  various  faces  of 
the  individual  sastrugi  (see  eq  13). 

The  resulting  model  does  quite  well  in  repro¬ 
ducing  the  values  of  CDNi0  that  we  observed  in 
AC95.  For  sastrugi  of  height  10  cm  and  with  an 
areal  coverage  of  15%,  the  model  predicts  that 
Cdnio  =  1*43  x  10-3  when  <X>  is  less  than  12°.  In 
AC95,  the  nominal  head-on  value  of  CDnio  that 
we  observed  was  1.5  x  10-3.  In  the  model,  as  O 
gets  larger  than  12°,  CDnio  increases  rapidly.  In 
this  example,  Cdnio  increases  by  17%  as  turns 
from  12  to  20°.  As  O  continues  increasing,  CDnio 
eventually  reaches  a  maximum  of  2.73  x  10"3  at 


O  =  101.3°.  This  increase  with  also  mirrors  the 
directional  sensitivity  we  reported  in  AC95, 
though  the  predicted  maximum  is  about  7% 
higher  than  our  measured  maximum  for  Cdnio/ 
2.54  x  10-3. 

The  upshot  of  this  modeling  is  that  the  small, 
10-cm-high  roughness  elements  associated  with 
drifting  snow — rather  than  pressure  ridges  or 
other  mesoscale  features — determine  the  local 
drag  coefficient  over  snow-covered  sea  ice. 
Though  some  have  stated  this  conclusion  explic¬ 
itly  (e.g.,  Joffre  1982)  or  implicitly  (e.g.,  Banke  et 
al.  1976, 1980),  my  model  is  the  first  to  provide  a 
theoretical  foundation  for  it. 

These  model  predictions  and  the  observations 
in  AC95  suggest  how  we  should  attempt  to  para¬ 
meterize  Cdnio  ha  terms  of  what  we  can  measure. 
Because  Cdnio  depends  crucially  on  how  the 
mean  wind  is  oriented  with  respect  to  the  sas¬ 
trugi,  the  key  is  predicting  sastrugi  orientation. 
Our  observations  showed  that  this  orientation 
depends  on  the  history  of  the  surface  wind.  By 
inferring  surface  winds  from  analyzed  fields 
of  the  geostrophic  wind,  we  should  be  able  to  de¬ 
duce  this  orientation.  The  event  to  look  for  is  a 
directionally  constant  surface  wind  of  at  least 
8  m/s  that  persists  for  12  hours.  At  the  end  of 
such  an  event,  Cdnio  will  be  roughly  1.5  x  10~3. 
Figure  4  then  shows  how  to  estimate  Cdnio  when 
the  wind  direction  subsequently  changes. 
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